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Abstract 



We discuss neutrino masses from higher than d — 5 effective operators in a supersymmetric 
framework, where we explicitly demonstrate which operators could be the leading contribu- 
tion to neutrino mass in the MSSM and NMSSM. As an example, we focus on the d = 7 
operator LLH u H u H d H u , for which we systematically derive all tree-level decompositions. 
We argue that many of these lead to a linear or inverse see-saw scenario with two extra neu- 
tral fermions, where the lepton number violating term is naturally suppressed by a heavy 
mass scale when the extra mediators are integrated out. We choose one example, for which 
we discuss possible implementations of the neutrino flavor structure. In addition, we show 
that the heavy mediators, in this case SU(2) doublet fermions, may indeed be observable at 
the LHC, since they can be produced by Drell-Yan processes and lead to displaced vertices 
when they decay. However, the direct observation of lepton number violating processes is 
on the edge at LHC. 
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1 Introduction 



From neutrino oscillations, it is evident that neutrinos are massive, see, e.g., Refs. [1,2]. 
If the neutrino masses originate from physics beyond the Standard Model (SM) and are 
suppressed by a high energy scale, it is convenient to parameterize the impact of the heavy 
fields, present in the high-energy theory, by a tower of effective operators O d of dimension 
d > 4. These operators made out of the SM fields, are invariant under the SM gauge 
group [3,4] (see also Ref. [5]). The operator coefficients are weighted by inverse powers of 
the scale of new physics A NP : 

£ = £ S M + ^f + £ c t 6 + "--, with £f ff0 c-^(9 d . (1) 

A NP 

Some of these effective operators result in corrections to the low-energy SM parameters and 
in exotic couplings. It is also known that there is only one possible operator at the lowest 
order in the expansion, £fg 5 , namely, the famous Weinberg operator [3], 

O w = (L c it 2 H) (Hir 2 L) (2) 

which leads, after Electroweak Symmetry Breaking (EWSB), to Majorana masses for the 
neutrinos. Here L and H stand for the SM lepton and Higgs doublets, respectively. At tree 
level, Ow can only be mediated in three ways [6] : by a singlet fermion, a triplet scalar, or a 
triplet fermion, leading to the famous type I [7-10], type II [11-16], and type III [17] see-saw 
mechanisms, respectively (see also Ref. [18]). For recent discussions of SUSY versions see for 
example [19-21] and references therein. Compared to the electroweak scale, the mass of the 
neutrinos in all three cases appears suppressed by a factor t>/A NP , where v is the Vacuum 
Expectation Value (VEV) of the Higgs boson. Substituting typical values, one obtains that 
the original see-saw mechanisms are pointing towards the Grand Unified Theory (GUT) 
scale. 

More recently, however, scenarios in which Anp ~ TeV have been drawing some atten- 
tion, since they are potentially testable at the LHC. In these cases, additional suppression 
mechanisms for the neutrino masses are required, and several possibilities open up: For 
example, the neutrino mass may be generated radiatively, where the additional suppression 
comes from loop integrals, or the smallness of the neutrino mass may be protected by lep- 
ton number, such as in the inverse see-saw. In this study, we instead argue that the d — 5 
operator in Eq. (1) is forbidden, and neutrino masses originate from higher dimensional 
operators [22-33] (see also Refs. [34-39] for related discussions). Note that there may be 
additional suppression mechanisms at work in this approach, such as loop suppression or a 
small lepton number violating parameter, see Ref. [28] for an example. 

There are several key ingredients to neutrino masses from higher dimensional operators [28] . 
Consider, for instance, the operator 

O 7 = {LLHH){rfH), (3) 

where we have omitted spin, flavor, and gauge indices. In this case, the (H^H) component 
can be closed in a loop, which means that the d — 5 operator is generated radiatively, and 
the d = 7 operator may not be the leading contribution to neutrino mass (depending on the 
new physics scale). We do not consider such operators in this work, which means that [28] 
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1. We have to forbid the lower dimensional operators by a U(l) or discrete symmetry. 

2. We need new (scalar) fields to construct the higher dimensional operators, since (H^H) 
is a singlet under any such symmetry. 

The simplest possibilities to enhance the field content of the SM are the addition of a Higgs 
singlet [23,24] 

C d = n+5 = -r^-A LLHH )(SY\ n = 1,2,3,... (4) 

A NP 

or the addition of a Higgs doublet, leading to the Two Higgs Doublet Model (THDM) [22, 
25,28,40] 

£cff 2n+5 = -^- 4 (LLH u H u )(H d H u ) n , n = 1,2,3,.... (5) 

A NP 

As it has been demonstrated in Ref. [28] in the framework of the THDM, the decomposition 
of Eq. (5) often leads to a linear or inverse see-saw structure [41-43] if two extra fermion 
singlet fields Nr and N' L are involved. The neutral fermion mass matrix then reads, in the 
basis {yl N R N' L C ), 

( (Yj)v e(Y?)\ 
M v = {Y v )v n> A NP . (6) 
Wl) A^ p fj," J 

Here e, //, and /i" are typically introduced ad hoc as small parameters, because they are 
protected by lepton number (for the e-term, see also Refs. [18,44]). However, if the neutrino 
mass is generated by a higher than d — 5 effective operator, some terms in the neutral 
fermion mass matrix in Eq. (6) can only originate from non-renormalizable interactions, 
which means that extra fields are needed and that the lepton number violating parameters 
are suppressed by powers of Anp- A general discussion and two specific examples leading 
to the pi- and e-term in Eq. (6) can be found in Sec. 3 of Ref. [28]. 

In principle, all neutrino mass models can be supersymmetrized. In practice, however, it 
turns out that supersymmetry often requires additional particles for consistency, mainly 
due to the boundary condition that the superpotential has to be holomorphic. In this work, 
we apply Ref. [28] to the framework of SUSY, more specifically, we start with the minimal 
supersymmetric standard model (MSSM) and the next to minimal one (NMSSM). A related 
discussion in the NMSSM framework can be found in Ref. [24]. In Sec. 2, we systematically 
discuss higher dimensional effective operators including two Higgs doublets and one scalar 
within SUSY. Then in Sec. 3, we focus on one specific decomposition leading to a linear or 
inverse see-saw in the form of Eq. (6). We will show that this requires an extension of the 
particle content of the model, which is potentially observable at the LHC. The main features 
will be outlined using a specific model, where we also discuss constraints due to existing 
data. In Sec. 4, we show that lepton flavor mixing related to neutrino physics is potentially 
testable at the LHC. We also show that the cross sections for some lepton number violating 
processes can be significantly larger than naively expected. However, it turns out that they 
are on the edge of observability at the LHC. Finally, in Sec. 5 we draw our conclusions. In 
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Op.# 


Effective interaction 


Charge 


d = 


5 


1 


LLH U H U 


2<ZL + 2q Hu 


d = 


7 


3 


LLH u H u H d H u 


2q L + 3qH u + QH d 


d = 


9 


7 


LLH u H u H d H u H d H u 


2q L + 4q Hu + 2q Hd 



Table 1: Effective operators generating neutrino mass in the MSSM up to d = 9. The operator numbers 
have been chosen in consistency with Table 2. 







Op.# 


Effective interaction 


Charge 


Same as 


d = 


5 


1 


LLH U H U 


2q L + 2q Hu 




d = 


6 


2 


LLH U H U S 


2ql + qu u - qu d 




d = 


7 


3 


LLH u H u HdH u 


2q L + 3q Hu + qu d 








4 


LLH U H U S S 


2qL - 2q Hd 




d = 


8 


5 


LLH u H u HdH u S 


2QL + 2q Hu 


#1 






6 


LLH U H U S S S 


2QL + 2q Hu 


#1 


d = 


9 


7 


LLH u H u HdH u HdH u 


2q L + Aq Hu + 2q Hd 








8 


LLH u H u H(iH u SS 


2qL + qH u ~ qu d 


#2 






9 


LLH U H U S S S S 


2<?l + qu u - qH d 


#2 



Table 2: Effective operators generating neutrino mass in the NMSSM up to d — 9. Here S is the NMSSM 
scalar, which means that its charge qs is fixed by the terms \SH u Hd, i.e., qs = —(qH u + QH d ), and kS 3 , 
i.e., 3qs = 0, which have been used to derive the charge condition (c/., Eq. (7)). 

Appendix A we comment on possible fundamental theories leading to the d = 7 operator 
discussed in this paper. In Appendix B we approximate some of the couplings for the model 
considered in Sec. 4. 

2 Neutrino masses from higher than d=5 operators in SUSY 

We take the MSSM as a starting point for various extensions. Note that its Higgs sector 
has the structure of a type II-THDM with the restriction that it is CP-invariant at leading 
order. As in Ref. [28], we require a discrete symmetry (in the sense of a matter parity) 
to forbid the d = 5 operator as leading contribution to neutrino mass, where the simplest 
possibility is Z 3 in the case of SUSY. 1 In Table 1, we list all possible higher dimensional 
operators made from lepton doublets and the two Higgs fields up to dimension nine. Note 
that compared to the THDM, the holomorphicity of the superpotential implies that the 
only possible effective d = 7 operator in the MSSM is LLH u H u H d H u , and it also limits the 

1 Note, that the requirement of forbidding the d = 5 operator automatically implies conserved R-parity 
because a AL = 1 operator immediately implies a contribution to the d = 5 operator [45]. For example, the 
sneutrino can get a VEV if R-parity is broken. Because of the neutrino-sneutrino-neutralino interaction, an 
additional d — 5 effective operator, which contributes to the neutrino mass, is then possible [46]. Therefore, 
we require R-parity conservation, which also has strong constraints on the possible decompositions of the 
effective operators. 
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number of possible decompositions further. 2 In Table 1, we show in addition to the possible 
operators the required charge combination such that the corresponding operator respects 
the discrete symmetry. However, in the MSSM, the combination fiH u H d appears in the 
superpotential, which is breaking such the discrete symmetry explicitly. For this reason, we 
consider models with extended particle content. 

A possible extension beyond the MSSM is the NMSSM, where an additional Higgs sin- 
glet S is introduced, see e.g. Ref. [48,49] for reviews. This singlet couples to the usual 
Higgs doublets H u and H d and obtains a non-zero vacuum expectation value. The NMSSM 
superpotential is 

WWssm = W Ynk + \SH u H d + kS 3 , (7) 

where HVuk denotes the superpotential for the Yukawa couplings, i.e., the MSSM superpo- 
tential. In this case, the charge of qs is fixed by the second and third terms being uncharged. 
From the discrete symmetry point of view, one can easily see that the last term in Eq. (7) 
is invariant under the Z3 for any charge assignment. 

In Table 2, we list all possible higher dimensional neutrino mass operators made from 
lepton doublets and the two extra Higgs fields up to dimension nine for the NMSSM. In 
the column "Charge" we also show the discrete symmetry charge using the fact that the 
terms in Eq. (7) have to be uncharged. In the last column "Same as", we indicate if the 
same condition as for a lower dimensional operator is obtained, i.e., the lower dimensional 
operator cannot be avoided in this case. One can read off the table that operators #2, #3, 
#4, and #7 can be independently chosen as leading contribution of neutrino mass, while 
the lower dimensional operators are forbidden. In general one can show that the NMSSM 
operators of the type LLH u H u (H d H u ) n S k (with n > 1, k > 1 or n — 0, k > 3) always 
imply that other operators of lower dimension are allowed as well. This is due to the fact 
that one finds field products of the type H u H d S or S* 3 , which have to be singlets under the 
discrete symmetry, since they appear in the superpotential Eq. (7). This means that d > 7 
effective operators generating neutrino mass with singlet scalars will always come together 
with lower dimensional operators. On the other hand, the effective operators with lepton 
and Higgs doublets only (such as #3 and #7) are per se interesting alternatives because 
one can choose even higher dimensional operators d > 9 as leading contribution. 

Note that operators #1, #2, and #4 have been studied in Ref. [24], whereas we focus on 
the d = 7 operator #3 in the following. In this case, a possible charge assignment for the 
Z 3 symmetry is 

<1h u = 0, q Hd = 1, q L = 1, (qs = 2) ■ (8) 

While in this case, both the MSSM and NMSSM can be used as a framework, one has to 
be aware of the fact that the /i-term of the MSSM, fiH u H d , explicitly breaks the discrete 

2 The holomorphicity of the superpotential implies that interactions among scalars of the form 4>i4>j'Pk 
can only be introduced via F-terms (another possibility to get non-holomorphic terms for the neutrino mass 
operator are non-canonical terms in the Kahler potential, as discussed in [37,47]). Since there are no SUSY 
invariant interactions with both, fermions and F fields, the only possible effective d = 7 operator in the 
MSSM is LLH u H u H d H u . 
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symmetry. This problem is automatically circumvented in the NMSSM, as Eq. (7) respects 
the discrete Z 3 symmetry and generates the //-term when S takes a VEV. In SUSY, there 
are some differences compared to the THDM case in Ref. [28]. For instance, the Lagrangian 
in Ref. [28] was invariant under a new U(l) symmetry in some cases, taking the role of the 
Z 3 symmetry here, which potentially lead to unwanted Goldstone bosons; see discussion in 
Sec. 3.1 of Ref. [28]. Even if the Lagrangian was invariant under such a symmetry, it is 
obvious that Eq. (7) would break it explicitly, while it respects Z 3 . In addition, note that 
a d — 5 operator is inevitably generated by connecting the external H d and H u lines of 
a d — 7 operator using a discrete symmetry breaking term m 2 Hd • H u (see discussion in 
Sec. 3.1 of Ref. [28]). The term fiH u ■ Hd in the superpotential corresponds to the scalar 
terms \fj,\ 2 H^H u and \/i\ 2 H d Hd in the Lagrangian here (see e.g. Sec. 16 of Ref. [50]), which 
means that this problem does not occur. Instead, one has a SUSY soft breaking term 
BfiH u ■ H d , which however would break the discrete symmetry and thus occurs only below 
Anp- Therefore it should be sufficiently smaller than Anp resulting in a suppressed d — 5 
one-loop contribution. Note however, that its value is bounded from below due to searches 
for the MSSM Higgs boson as it is proportional to the mass of the pseudoscalar Higgs boson. 

The possible decompositions of LLH u H u HdH u are systematically derived in Appendix A at 
tree level, where the mediators for different possibilities are listed in Table 6. Note that the 
right-handed fields listed there have to be incorporated as charge-conjugated left-handed 
fields in the superpotential. These decompositions can be roughly categorized as extensions 
of the well known d — 5 decompositions, the type I, II, and III see-saw scenarios. We define a 
decomposition as extended type II see-saw, if all mediators are scalars, i.e., decompositions 
#5, #6, and #21-^24 in Table 6. Therefore the only appearing fermions are the external 
lepton doublets. The only lepton number violating interaction is then 

(77ir 2 fL)0, (9) 

where represents one of the scalar mediators. This vertex violates lepton number by 
AL = 2 and therefore conserves R-Parity. All other decompositions that have fermionic 
mediators can be seen as extensions of the d = 5 type I or type III see-saw mechanism. 
Since we can have several combinations of scalar fields and SU(2) singlet, doublet, or triplet 
fields as mediators, a further distinction will not be made. Depending on the topology 
and the actual realization of these operators, the various decompositions have different 
characteristics. Note that integrating out all but two neutral fermion fields will lead to an 
inverse see-saw-like scenario, as in Eq. (6). 

Let us now illustrate some of the complications involving extra scalars as mediators. As an 
example we take decomposition #1 from Table 6, shown in Fig. 1. One can easily see that 
the scalar <fi has the same quantum numbers as the scalar of the NMSSM, and it also has the 
same coupling to the Higgs fields. If the terms in Eq. (7) are present in the superpotential, 
it can get a VEV iy This in turn means that the d = 6 operator of the type LLH U H U S (#2 
in Table 2) may be allowed, where <f> = S, which may be the leading contribution to neutrino 
mass. Indeed one can see from the XHuHdft -vertex in Fig. 1 that we have for the discrete 
symmetry charge = qn u + QH d , which means that we cannot forbid the d = 6 operator 
LLH U H U S which leads to neutrino mass if <fi obtains a VEV. In summary, the MSSM 
extended by a scalar singlet mediator can potentially be NMSSM-like, and can potentially 
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Figure 1: Decomposition #1 from Table 6 of the effective d = 7 operator LLH u H u HdH u . Here N and 
N' are fermion singlets, and is a scalar singlet. 




Figure 2: Decomposition #17 from Table 6 of the effective d = 7 operator LLH u H u HdH u . Here iV and 
N' are fermion singlets, and £ and £' are fermion doublets. 

induce the d — 6 operator, which may dominate the neutrino mass contribution. Since the 
operator LLH U H U S in the NMSSM has been studied in Ref. [24], and since we want to avoid 
the d — 6 operator genuinely, we focus on decompositions with two neutral fermions (to 
reproduce the inverse see-saw) and no singlet scalars in the following. One of the simplest 
examples is decomposition #17 in Table 6, which we will discuss in greater detail, see 
Fig. 2. While for neutral SM singlets as mediators the production rates of the new particles 
are rather low, the SU(2) doublets in #17 will lead to gauge interactions with potentially 
observable phenomenology at the LHC. However, note that also the fermion singlets could 
be replaced by triplets, which would lead to a see-saw Ill-type phenomenology. 



3 A linear or inverse see- saw example 



We have seen in the previous section that extensions of the MSSM containing NMSSM- 
like singlets do have some problems. Therefore we consider a model where we add two 
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gauge singlet superfields N and N' and a vector-like pair of SU (2) doublets £ and £' with 
hypercharges Y — \ and F = — |, respectively. Note, that the holomorphicity of the 
superpotential requires to use left-handed fields only and, thus, £ corresponds to the charge 
conjugated field 2 R l ^ 2 of decomposition #17 in table 6. Moreover, none of these additional 
fields can participate in the breaking of the electroweak symmetry as we require R-parity 
conservation. The corresponding superpotential is 

W = W quavh} + Y e e c L ■ H d - Y N NL ■ H u + Kl N'£ ■ H d - k 2 N'£' ■ H u + m N NN' 

, . . . (10) 
+ ■ £ + ^H u ■ H d , 

and the corresponding lepton number assignments are 

L(N) = -1 , L(N') = +1 , L(i) = -1 , L(i') = +1 (11) 

implying that the interaction proportional to k 2 breaks lepton number by two units 3 . This 
superpotential yields the following part of the Lagrangian for the fermions carrying lepton 
number 

£ fcrmionic = _ . ^ + ^ L . ^ + ^ . fi d ) + Y N (NL • H u + NL • H u + NLH U ) 

- Kl (N'S ■ H d + N't ■ H d + N'l • H d ) + k 2 (N'C' ■ H u + N't? ■ H u + Nf ■ H u ) 

- m N N'N - m£' ■ £ + h.c. , 

(12) 

using 2-component spinors. For completeness, note that for the leptons we use the "~" for 
scalars, whereas for the Higgs bosons we use the "~" for the fermionic partners. After the 
Higgs fields get a VEV, the mass matrix for the neutral fermions reads in the basis 

f° = (u,N,N',e,e°) (13) 



/ Y N v u \ 

F>„ m T N 

m N Kiv d k 2 v u 

njv d — 

\ kJvu —m^ J 



(14) 



The corresponding mass eigenstates will be denoted by rii with \rrii\ < \m,j\ for i < j. The 
mass terms for the charged fermions are given by 

-v d e c Y e e L - m^ + £'- . (15) 

We can now determine the neutrino mass, by integrating out the mediator fields. For the 
sake of simplicity, let us first of all ignore the flavor structure. Integrating out the heavy 
doublets, we obtain 

C N = N^a^d.N + N'ha^d.N' - m N NN' - Y N NL ■ H u - — N'N'H U ■ H d + h.c. , (16) 



3 Note that this assignment is, to some extent, arbitrary, and that a different assignment couid be chosen 
such that the interaction proportional to K\ breaks lepton number by two units. However, none of our 
conclusions is affected by this specific choice. 
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which reads in the basis {y, N, N') after electroweak symmetry breaking 






Y N v u 





Y N v u 





m N 





m N 





M° f = ( Y N v u m N | (17) 

with jl = v u v d {2K 1 K 2 )/m^. This is an inverse see-saw mass matrix, as in Eq. (6). As 
it is characteristic for the inverse see-saws from higher dimensional operators, the lepton 
number violating term is suppressed by a heavy scale. If in addition the singlet fermions 
are integrated out, we obtain for the neutrino mass scale 

m v = vlv d Y*^ r . (18) 
mgmr N 

For a neutrino mass m v ps 1 eV and v ~ 177 GeV and the heavy mass scale at ITeV this 
means couplings O(10~ 3 ) are required. Note that this coupling strength in not unreasonably 
small, although these couplings are, in addition, protected by lepton number. 

If we assume instead a hierarchy of the heavy particles where the isospin singlets are heavier 
than the doublets, we first can integrate out the singlets. The mass matrix for the remaining 
neutral fields (u, reads then 



Mf = \R lVd -m 5 | , (19) 
K k 2 v u — m € 

where kx/ 2 = «i/2 Y^/m N . Integrating out the £ fields afterwards we again arrive at Eq. (18) 
for the mass of the light neutrinos. 

As in the conventional inverse see-saw, in which the /t-term in Eq. (17) is introduced ad hoc, 
there are various interesting phenomenological effects of this scenario. It is expected that 
non-unitarity and its CP violation can be tested at possible long-baseline neutrino oscillation 
experiments (see, e.g., Refs. [51,52]). Furthermore, one may observe lepton-flavor-violating 
(LFV) processes such as fj, — > ey. Lepton-number-violation, on the other hand, is expected 
to be hardly testable in conventional scenarios, since the heavy Majorana neutrinos form 
pseudo-Dirac particles with suppressed Majorana character, see, e.g., Ref. [53] 4 . 

There are several ways to realize a flavor structure that is in accordance with neutrino 
physics. Since there are three distinct (active) mass eigenstates, at least two of them must 
be massive. The straightforward approach is to add three generations of the heavy fields, 
which leaves, however, many unconstrained parameters. Another possibility is to generate 
one neutrino mass by the inverse seesaw with one generation of mediators, and the second 
neutrino mass at the one-loop level [57] if the flavor structures in the soft SUSY sector 
differs from the ones in the superpotential. A third version is the minimal inverse seesaw 
scenario in Ref. [58], consisting of only two generations of the heavy fields, which narrows 
down the number of free parameters. 

We follow a similar approach, where we assume that one neutrino state is massless. We 
assume two generations of N and N' each, and only one generation for the other mediators. 



ki v d k 2 Vu 



'Some attempts to avoid the suppression are discussed in e.g., Refs. [54-56]. 
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Thus, compared to Eq. (18), we obtain a mass matrix 

(rn v ) aP = vlv d (Y N ) ai (mJ f 1 ) ij /j, jk (m N 1J ) k i(Y^)i l3 , 

where 

1 



Pjk 



((«l)j(K 2 )fc + (« 2 )j(«l)fc) 



(20) 



(21) 



The flavor basis can be chosen in a way that M N (and consequently M^ 1 ) is diagonal, with- 
out loss of generality. We choose the parameters to reproduce tri-bimaximal mixings [59] 5 , 
which depend on the mass hierarchy: 

Normal hierarchy. A rather straightforward parameterization is: 



Y N = y N 



\V3 



\ 

V2 ) 



K\ = ki 



-1 
1 



k 2 = k 2 



where 



P = A/m 2 /m 3 , 2vlv d y 2 N k 1 k 2 /(m 2 N m^) = m 2 . 



(22) 



(23) 



This reproduces the tri-bimaximal mixing pattern and two non-zero mass eigenvalues. In 
this case, the flavor structure of the neutrinos is dominantly generated by the neutrino 
Yukawa couplings Y N . Since we have more parameters than constraints from neutrino 
physics, there is a certain freedom in the parameters of the couplings. For example one can 
vary y^ as long as this is compensated by an according change of or k\/ 2 . The mass 
ratio p = \fm 2 /m- i can be generated by y^, or k\i 2 . A possible set of parameters is 
3y N = 10" 3 , ki = k 2 = 10" 2 , tan/3 = 10, m N = 1070 GeV, and m 5 = 200 GeV, which we 
will use in the next section. We have checked that this point is compatible with bounds 
on rare lepton decays such as p — > cy as well as with the search for the tri-lepton signal of 
supersymmetric particles at the Tevatron [61] and searches for new physics in final states 
containing leptons at the LHC [62,63]. In order to satisfy the bounds from the rare decays, 
we have assumed that the scalar leptons are so heavy that their contributions are suppressed 
and that the leading contributions are due to loops containing fermions and the W-boson. 
Note that the product v^Vd in Eq. (23) peaks at about tan/3 ~ 2, and it becomes small for 
large tan/3. 

Inverted hierarchy. The inverted hierarchy can be obtained by the parameterization 



Y, 



N 



Vn 



( 



2 
1 
1 



f , 

/3/ 



Ki = ki 



k 2 = k 2 



M N = m N 



(24) 



where p = \/mi/m 2 and 



2vlv d y%kik 2 / (m^rrii:) = m 1 



(25) 



In the following, we only consider the normal hierarchical example. 



3 If #i3 > 0, as indicated by the recent T2K hint [60], a different flavor structure can be easily implemented. 
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Figure 3: Total cross section a(pp — > ^f ) as a function of the mass m^. 

4 LHC phenomenology 

In many supersymmetric versions of neutrino mass models one finds traces of the underlying 
mechanism generating neutrino masses in the spectrum and decay properties of the super- 
symmetric particles, for an incomplete list see e.g. [19,21,64-74]. Before discussing this in 
more detail, let us have a look on the number of parameters related to neutrino physics in 
this model. Working a basis where the lepton Yukawa couplings are flavor diagonal, the 
following parameters contribute: Yjv, and which amounts in our specific model 

into 24 real parameters if all CP violating phases are taken into account, but which gets 
reduced to 13 if all phases are zero. From these at most six can be determined in the near 
future, leaving 18 (7) parameters undetermined. Here the question arises to which extent 
they might be measured or at least constrained at the LHC. In principle we have sufficient 
many decays to determine them: six heavy neutral states decaying into the three charged 
leptons, in total 18 decays. 

Of particular interest is the question to which extent the new particles, which are postulated 
in Eq. (10), can be produced at the LHC. It turns out that except for finely tuned parameter 
combinations, the heavy states are either mainly SU(2) singlets or mainly SU (2) doublets 
corresponding to the electroweak states. The singlet states contained in the superfields N 
and N' can only be produced in cascade decays. However this will happen in rare occasions 
only due to the smallness of the involved Yukawa couplings. The SU (2) L doublets contained 
in £ and on the other hand, can in principle be produced directly in Drell-Yan processes 
similarly to sleptons or charginos and neutralinos within the usual MSSM [75]. In Figure 3 
we show the cross section for as a function of for y/s = 7 TeV and 14 TeV assuming 
that the mixing with the singlet fields is small. The numbers have been obtained using 
the WHIZARD package [76]. If one takes for example a mass of 200 GeV for £ + and £° 
using the numerical values for the couplings as given in the previous section we find for 
the total cross section cr(pp — > £ ± £°) about 122 fb (417 fb) in case of 7 (14) TeV cms- 
energy using the WHIZARD program [76]. Note, that the states and are mainly a nearly 



10 



maximal mixed superposition of the neutral components of the SU (2) doublets in this case 
n 4 , 5 ^(£°±0/v/2. 

In the following, we concentrate on the fermionic states, as they are directly related to 
neutrino physics. We assume for the moment that their scalar partners are much heavier, 
so that the only possible decay channels are into SM fermions and Higgs bosons, or vector 
bosons. In this case, the main decay modes are 

e -> W + v k , H + u k (26) 

for £ + , which is the Dirac fermion composed of the charged components of £ and Note, 
that there are no decays into Z or h° as this particle does not mix with the charged leptons 
at tree-level. One expects small decay widths as these decays have their origin in the mixing 
of £° with neutrinos and, thus, their widths are proportional to the corresponding mixing 
matrix element squared. We indeed find = 1.42 • 1CT 5 keV. In case of the neutral 

fermions rii, a larger variety of decay channels is possible: 

m -> W ± IJ, HHJ (27) 
7n ->■ Zu k , h°v k , H°v k , A°v k (28) 

with lj = e, (/,, t. In these cases, the decays also originate from the mixing of the neutral 
states with the neutrinos and, thus, the corresponding widths are expected to be small as 
can be seen in Table 3 where we give the corresponding widths and branching ratios for the 
scenario discussed in the previous section. Obviously some of the widths are so small that 
one can expect sizable decay lengths at the LHC in the range of 100 /im to several mm once 
the boost factor is taken into account. This is an important feature because in this way 
one can not only suppress the SM background, but one can also identify the leptons coming 
from these decays and distinguish them from the ones coming from the cascade decays of 
supersymmetric particles. Another interesting feature is, that there are two pairs of states 
where within each pair the branching ratios are nearly equal: n^jn-j and n s /ng. The reason 
is that they form a quasi Dirac fermion. Also in case of n 4 and n 5 we have a quasi Dirac 
fermion resulting in difficulties to determine the branching ratios of the individual states. 
As a consequence at most 9 branching ratios can be related to neutrino physics in praxis. 

The fact that the rii decay into W ± l T clearly proves that these particles carry lepton number 
and, thus, one might suspect that they are related to the generation of neutrino masses. An 
important question is in this context to which extent one can prove their Majorana nature 
by observing both lepton charges in the final states. Therefore one has to look for final 
states which violate lepton number by two units compared to the initial one: 



ud 


-> 


z+r 


~w~ 


(29) 


ud 


->■ 


i + t 


-w~z 


(30) 


qq 


->■ 


i + t 


-w~w-,ri'~w + w + 


(31) 



Note that the leptons can easily be of different generations due to the large mixing angles 
in the neutrino sector. In the calculation of these processes we have included all possible 
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Particle 


T[keV] 


BR{W ± e T ) 


BR(W ± fi T ) 


BR(W ± r T ) 


BR(Zv) 


BR(h°u) 




2.3 • KT 5 


6.6 • 1(T 3 


7.0-10- 2 


0.18 


0.36 


0.38 


n 5 


1.9 • 10 -5 


1.2 • 10 -2 


0.41 • 10- 2 


0.18 


0.42 


0.34 


n & 


1.2 


1.2 • lir 11 


0.21 


0.21 


0.21 


0.37 


n 7 


1.2 


1.2 • KT 11 


0.21 


0.21 


0.21 


0.37 


n 8 


2.9 


0.14 


0.14 


0.14 


0.20 


0.38 


n 9 


2.9 


0.14 


0.14 


0.14 


0.20 


0.38 



Table 3: Total decay widths of the neutral mass eigenstates and branching ratios into the possible final 
states (where v is the sum over the three light neutrino mass eigenstates). 



intermediate particles to account for possible mixing effects, e.g., due to the pseudo-Dirac 
nature of the heavy additional neutral fermions, which turn out to be important in case of 
lepton number violating final states. The interesting part of the Lagrangian is given by 

£wi jni = -^n a ^a l3 P L + k 3 P R )l 3 W+ - JLl^a*^ + b^P^W' (32) 
Vn, = —jf t Y{c t P L + diP R )CW+ - ^=tr(c*P L + d*P R ) nj W- (33) 



where 



Oij = Uij , bij = j = e,fx,r, (34) 



since these couplings have their origin in the left-handed couplings of the SM leptons to 
the J¥-boson. Here we are working in a basis where the Yukawa matrix of the charged 
leptons is diagonal. Note, that would only be non-zero if there were a mixing between 
£- with the charged leptons. The couplings to the £~ = (£ + )r) T originate from the 
SU (2) couplings between the £ doublet components and are given by 

d = Ui? , di = -U*£ . (35) 

Here U denotes the matrix diagonalizing the mass matrix of the neutral fermions, see 
Eq. (14). 

The results for the 2 — > 3 processes are shown in Table 4. Note that in this case only final 
states containing a H^-boson are possible, as the £ + does not decay into charged leptons. 
The main contributions in this case are due to 

ud ->■ l + n* (i = 1, . . . 5) (36) 

as shown in Fig. 4. Here the rij are in this case either mainly neutrinos (i = 1,2,3), or 
an admixture of £° or £ 0/ (i = 4,5). The contributions of the neutrino-like states are 
suppressed because they are off-shell, whereas in case of the £°/£°'-like states, there are on- 
shell contributions. These are, however, somewhat suppressed by the small mixing elements 
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Figure 4: Dominant contribution to the LNV processes ud — > W e + e + and qq — > e+e+. 



with the neutrinos. We have put a cut on the invariant mass of the leptons of 10 GeV as 
otherwise the e + e~ final states would be enhanced by several orders of magnitude due to an 
nearly on-shell photon. The flavor mixed final states are of the order of a few fb and thus 
are potentially observable if sufficient luminosity is accumulated. Note that for extracting 
the corresponding signal, only the hadronic final states of the VT^-boson should be considered 
and, thus, the cross section shown has to be multiplied by the corresponding branching ratio 
BR(W — > qq'). In the case that the two leptons have different flavor, these processes are 
essentially background free. However, in case of equal flavor leptons multi H^-production in 
association with a Z-boson or an off-shell photon will contribute. Both contributions can 
be suppressed by putting cuts on the invariant mass of the two leptons. 

One also sees from these tables that the lepton number violating final states are strongly 
suppressed which is due to the appearance of the pseudo-Dirac like state n^/n^ implying 
that the final contribution to the cross section is proportional to m^ 5 — m\ A ~ 0(ml), and 
thus tiny. 

In Table 5, we give cross sections for the 2 — > 4 processes containing two H^-bosons. Note 
that we do not give the corresponding ones containing a Z-boson, see Eq. (30), which 
are smaller because the corresponding contributions are those of 2 — > 3 processes plus an 
additional Z-boson, attached to all internal and external lines in case of the lepton flavor 
mixing/violating final states. As expected, the cross sections of lepton flavor conserving 
and lepton flavor mixing final states are about two orders of magnitude smaller than the 
ones of the corresponding 2 — > 3 processes. However, the cross sections for the lepton 
number violating processes are larger than naively expected. This can be understood as 
follows: in case of the 2 — 3 processes all lepton number violating contributions are due to 
the Majorana nature of the neutral fermions and are suppressed by the pseudo-Dirac like 
nature of the heavy states. In case of the 2 — > 4 processes, there are additional contributions 
which are proportional to the momentum of the off-shell neutral particles times two powers 
of lepton number violating couplings, e.g. they are proportional to |cjrfj| 2 . Performing an 
approximate diagonalization of the neutral fermion mass matrix as done in Appendix B, one 
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Process 


(j [fb] (7 TeV) 


a [fb] (14 TeV) 


nn V 

PP ^ 




"e 


(1.651 ± 0.024) • 10 2 


(4.161 ± 0.023) • 10 2 


Tin V 

PP ^ 






(9.240 ± 0.033) • 10 


(2.671 ± 0.042) • 10 2 


nn ^ 

pp ^ 


W + e + fi- 


(1.068 ±0.099) 


(2.848 ± 0.011) 


nn V 

PP ^ 


W + e~ 




(1.057 ±0.013) 


(2.871 ± 0.012) 


nn ^ 

pp ^ 


W~e + fi- 


(5.748 ±0.015) • 10" 1 


(1.742 ± 0.015) 


nn V 

PP ^ 


W~e~ 




(5.755 ± 0.015) • 10 _1 


(1.753 ± 0.017) 


nn ^ 

pp ^ 


W + e^ 




(1.058 ±0.096) 


(2.861 ± 0.011) 


nn ^ 

PP ^ 


W + e~ 


-T+ 


(1.056 ± 0.095) 


(2.854 ± 0.011) 


nn V 


W~e^ 


~T~ 


(5.714 ± 0.015) • 10 _1 


(1.754 ± 0.015) 


nn V 


W~e 


-T+ 


(5.750 ± 0.015) • 10 _1 


(1.744 ± 0.019) 


nn V 

PP ^ 




>~ 


(1.676 ± 0.014) • 10 2 


(4.116 ± 0.023) • 10 2 


PP ->■ 




>~ 


(9.242 ± 0.033) • 10 


(2.677 ±0.035) • 10 2 


PP ^ 




-T- 


(2 fifis + n oid) ■ i rr 1 

l^.\J\J(J _1_ \J . \J Zjt: 1 1U 


(7 092 ± 028^1 • 1 _1 


nn ^ 

pp ^ 




-T + 


(2.652 ±0.026) • 10 _1 


(7.187 ± 0.029) • 10" 1 


PP ^ 


W~fT 


-T- 


(1.432 ±0.006) • 10" 1 


(4.424 ± 0.038) • 10" 1 


PP ^ 


W~fi- 


"r+ 


(1.439 ±0.004) • 10" 1 


(4.433 ± 0.037) • 10" 1 


PP ^ 


W+T+T- 


(1.665 ±0.023) • 10 2 


(4.138 ±0.063) • 10 2 


PP ^ 




~r 


(9.265 ± 0.034) • 10 


(2.652 ± 0.035) • 10 2 


nn ^ 

pp ^ 


W~e^ 


"e+ 


(4.711 ± 0.069) • 10~ 12 


(4.847 ± 0.030) • 10~ n 


nn V 

PP ^ 


W + e 




(1.423 ± 0.008) • 10~ 12 


(1.818 ± 0.071) • 10~ n 


nn ^ 

pp ^ 


W~e^ 


> + 


(1.017 ±0.014) • 10~ n 


(9.869 ± 0.054) • 10~ n 


nn V 

PP ^ 


W + e~ 


/i" 


(3.184 ± 0.015) • 10~ 12 


(3.22 ± 0.15) • 10 -11 


nn V 

PP ^ 


W~e+T + 


(1.169 ± 0.015) • 10~ n 


(1.050 ± 0.054) • 10~ 10 


— > 


W + e~ 


"r 


(4.173 ± 0.020) • 10~ 12 


(4.12 ± 0.28) • 10~ n 


— > 


W-fi + /j + 


(5.861 ±0.082) • 10~ 9 


(2.278 ±0.013) • 10~ 8 


pp — >■ 


W + fi~ 


> 


(2.377 ±0.010) • 10~ 9 


(1.153 ±0.017) • HT 8 


pp — )■ 


W~fi + r + 


(1.184 ±0.013) • 10~ 8 


(4.584 ±0.023) • 10~ 8 


pp — )■ 




"r 


(4.788 ±0.018) • 10~ 9 


(2.363 ±0.039) • 10~ 8 


pp — >• 




-T+ 


(5.956 ±0.080) • 10~ 9 


(2.292 ±0.031) • 10~ 8 


pp — >• 


W + T- 




(2.383 ±0.010) • 10~ 9 


(1.120 ±0.014) • lO" 8 



Table 4: Cross-sections for the processes with W ± £ ± £ ± as final states (lepton number violating processes 
in lower section). A cut on the invariant lepton mass of 10 GeV has been assumed. 
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a [fb] (7 TeV) 


rr [fhl f14 TpVI 




pp — >■ 14/ + e 14/ 


(3.447 ± 0.87) 


■ 10 1 


(1.277 ± 0.66) 




pp —tW e w \i 


(7.06 ± 0.15) • 


10 3 


^o.l41 ± U.UZ I ) ■ 1U 


-2 


— > W e W li 


(6.99 ± 0.16) • 


10 13 


{o.Zvb ± U.Uz g • 1U 


-2 


\ 1/1/ — r~ ^ — 1/1/ — — r" 

pp — > W e W 


(1.037 ±0.020) 


• 10~ 2 


( A OOQ _J_ n HQfi^ 1 n- 

(^4.zyo ± U.UooJ • 1U 


-2 


TT7"4- -I-TT7" — — 

pp — > W + e + W t 


(1.015 ± 0.021) 


• 10~ 2 


(4.411 ± 0.036) • 10 


-2 


mm . TIT" + < T T * , + 

pp W^fi W ji 


(3.74 ± 0.10) • 


10 1 


(1.279 ± 0.017) 




Pj9 — > /i W T 


(2.913 ± 0.048) 


• 10 3 


(l.Uyo ± U.UU I ) ■ 1U 


-1 


pp — >■ /i r 


(2.990 ± 0.042) 


-i r\ — 9 

• 10 2 


/ 1 i on _i_ n nn7\ i n- 

(i.idy ± u.uu • iu 


-1 


pp — f 14 7 T W T 


(4.27 ± 0.10) ■ 


10 1 


(l.oUo ± U.U1 1 ) 




pp — >■ vv T e W^e 


/ 1 110 _i_ n n 1 o \ 

(1.112 ± 0.013) 


• 10 


(A o/^ 1 _i_ n noo\ 1 n- 

(4.Z61 ± U.Uz8) • 1U 


-4 


pp r Vvt^VV fX 


(1.537 ±0.023) 


• io- 3 


^ 8m + n n^n^ . i 

^o.oiu nz u.uouj iu 


-3 


pp ->■ W^+e'W^+r" 


(4.721 ±0.055) 


• 10~ 3 


(1.761 ±0.016) • 10" 


-2 


pp ->■ W + ltW + n~ 


(4.099 ±0.052) 


• 10~ 3 


(1.514 ±0.013) • 10" 


-2 


pp -+ W + Ll~W + T~ 


(2.704 ±0.036) 


• io- 2 


(1.062 ±0.093) • 10" 


-1 


pp ->■ W + T~W + T~ 


(4.614 ±0.065) 


• io- 2 


(1.729 ±0.016) • 10" 


-1 



Table 5: Cross-sections for the processes with W + £ W ± t f as final states (lepton number violating 
processes in lower section). A cut on the invariant lepton mass of 10 GeV has been assumed. 



sees that this combination of couplings does not vanish in the limit of vanishing neutrino 
masses as they are roughly proportional to 



(Qj/Ci + bjK 2 ) 

M 4 N mj 



(37) 



e.g. they only vanish in the limit where either one of the heavy masses goes to infinity or 
both couplings, Ki and k 2 , to zero. This is a consequence of the fact that £ and £' form a 
vector-like representation of SU(2). 

In summary we find that one should be able to detect the SU (2) doublets up to masses of 
about 1 TeV and show that they carry lepton number. However, it turns out that the cross 
sections for the processes violating total lepton number are on the edge to be discovered, 
as they would require at least a luminosity of the order of ab~ l in the most optimistic 
cases, e.g., by considering at least 10 events without any background considerations due to 
detector effects. 
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5 Summary and conclusions 



In this work we have studied neutrino mass generation from higher than d — 5 effective 
operators in supersymmetric models. While the d — 5 operator typically points towards 
the GUT scale, higher dimensional operators may be generated by mediators observable 
at the LHC If any d = 5 contribution is to be forbidden, a discrete symmetry is needed, 
which can be used to control the dimension of the effective operator generically dominating 
neutrino mass. While this discrete symmetry is to be softly broken by the //-term of the 
MSSM, the Z 3 symmetry, the SUSY Lagrangian is invariant under, can be naturally used 
in the NMSSM. We have also taken into account that in the NMSSM, higher than d = 5 
effective operators leading to neutrino mass may include the NMSSM scalar and the two 
Higgs doublets. While the NMSSM scalar can be used in d — 6 and d = 7 effective operators, 
for d > 7, only Higgs doublets are allowed in the effective operators, since otherwise lower 
dimensional effective operators are generated as well. Therefore, we have focused on the 
d = 7 operator LLH u H u HdH u as the simplest possible example in the following, which 
respects this line of argumentation. 

For this operator, we have derived the list of possible decompositions at tree level system- 
atically. The results have been similar to an earlier work [28] , with the exception that some 
topologies have been forbidden by the holomorphicity of the superpotential. Many of the 
derived decompositions can be regarded as extensions of the usual type I, II, or III see-saw 
mechanisms because of a similar field content. Models with two extra heavy fermion singlets, 
for example, lead to inverse see-saw scenarios if the additional mediators are integrated out, 
where the lepton number violating term is naturally suppressed by the mediator mass. As a 
peculiarity of supersymmetry, we have identified that singlet scalars are potentially harmful 
because they may induce lower dimensional operators dominating neutrino mass if similar to 
the NMSSM scalar. Therefore, we have chosen an example with two extra fermion singlets 
and heavy lepton doublets which are vector-like under SU(2). We have also demonstrated 
how the flavor structure for normal and inverted hierarchy can be easily implemented using 
two generations of the heavy fermion singlets. 

Focusing on the new fermions, we have demonstrated that parts of the model can already be 
tested with the current LHC run at 7 TeV by displaced vertices, and at 14 TeV we expect 
that it can be tested up to masses of several hundred GeV for the SU (2) doublets. We have 
also seen that the cross sections of some of the lepton number violating processes are larger 
than naively expected, but still on the edge of observability at the LHC. 
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Topology 1 Topology 2 Topology 3 Topology 4 



Figure 5: Possible topologies for the effective d — 7 operator LLH u H u HdH u . Topologies 3 and 4 cannot 
be realized in SUSY. Solid lines are either fermions or scalars, dashed lines are always scalars. 

A Possible decompositions for the operator LLH u H u HdH u 

In this appendix, we systematically discuss underlying, more fundamental models leading 
to operator ^3 in Table 1 or Table 2 at tree level. The results are similar to Ref. [28] in the 
THDM. The possible topologies for the decompositions can be found in Fig. 5. In SUSY, 
however, topologies 3 and 4 can be excluded. This is due to the fact that scalar couplings 
in SUSY have to be of the type cftcj), <ft<f><f>, or $ $ 4>4>i since they are generated by 

F-terms and D-terms, as a consequence of holomorphy. In topology 3, the scalar four-vertex 
has to be of the type HHX*X* , where X is a heavy virtual scalar field. The three- vertex 
must be HHX* . These two vertices can not be connected by a propagator A X - Hence 
topology 3 can not be realized in SUSY. In topology 4, the four-scalar vertex can only be of 
the type HdH u H u X^ or H U H U H U X^ in order to produce an effective operator of the type 
LLH u H u HdH u . The only possible scalar four couplings allowed by SUSY, however, are of 
the type 0^0^00. Hence also topology 4 is not possible. The most economical extensions of 
the (N)MSSM may use the superpartners of the SM fields as mediators. However, at least 
at tree level and with R-Parity conservation, this is not possible. As all external fields, L, 
H u and Ha, have R — +1, a mediator with R — — 1 would cause vertices where R-Parity 
is violated. As a consequence, we have to introduce additional fields as mediators, and we 
also obtain superpartners for them. The possible decompositions of the d = 7 operator 
LLH u H u H d H u at tree level are shown in Table 6, where the brackets refer to the vertices 
with external fields for any given topology. If r appears, the fields couple to a triplet 
mediator; if not, they couple to a singlet. The mediators are denoted by X y , where 

• X denotes the SU(2) nature, i.e., singlet 1, doublet 2, or triplet 3. 

• C refers to the Lorentz nature, i.e., scalar (s), vector (v), left-handed (L) or right- 
handed (R) chiral fermion. 

• Y refers to the hypercharge Y = Q — 1^ . 

Besides the fixed sign of the scalars' hypercharges and the forbidden topologies 3 and 4, the 
decompositions are similar to the THDM case in Ref. [28]. Note that R and L indicate right- 
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# Operator Top. Mediators 

1 {H u ir 2 T^){H u ir 2 L){H d ir 2 H u ) 2 1* 1 L , lg 

2 (HjT 2 TlS)(H u iT 2 L)(H d iT 2 TH u ) 2 3* 3 L 1« l£, 3g 

3 (H u iT 2 fL-){H u iT 2 fL)(H d ir 2 H u ) 2 3 R , 3£, lg 

4 (-ie a6c )(iJ u ir 2 r a L c )(F u ir 2 r fe L)(^ir 2 r c i7 u ) 2 3* 3£, 3g 

5 (L°ir 2 f L) (H d ir 2 H u ) (H u ir 2 rH u ) 2 3^, 3^, lg 

6 (-ie abc )(r c ir^ a L)(H d iT 2 T b H u )(H u ir 2 T c Hu) 2 3^, 3^, 3g 

7 (H u ir 2 L c )(LiT 2 rH d )(H u ir 2 fH u ) 2 1* lg-, 3 R X , 3£ l5 3^ 

8 (-ie ab %H u irYL c )(Lir 2 r b H d )(H u ir 2 T c H u ) 2 3£, 3#, 3^, 3^, 3 S +1 

9 (^ u ir 2 L-)(ir 2 /J u )(L)(/J d ir 2 if u ) 1 1* l£, 2* /2 , 2* 1/2 , lg 

10 (H u ir 2 TL-)(ir 2 rH u )(L)(H d iT 2 H u ) 1 3ff, 3j, 2* /2 , 2^, lg 

11 (H u \T 2 ~Lf) {it 2 H u ){fL) {H d iT 2 fH u ) 1 1 R , l£, 2* 1/2 , 2^ 1/2 , 3g 

12 (H u ir 2 r a L-)(iT 2 r a H u )(T b L)(H d iT 2 T b H u ) 1 3 R , 3£, 2 fi 1/2 , 2 L 1/2 , 3g 

13 {H u ir 2 lJ){L){iT 2 H u ){H d iT 2 H u ) 1 1* l£, 2 s +1/2 , lg 

14 (/J u ir 2 rL c )(fL)(ir 2 iJ u )(i7 (i ir 2 /J u ) 1 3* 3£, 2 s +1/2 , lg 

15 (H u ir 2 L c )(L)(ir 2 rH u )(H d ir 2 fH u ) 1 1 R , l£, 2« +1/2 , 3g 

16 (H u ir 2 r a T-)(j a L)(ir 2 r b H u )(H d ir 2 r b H u ) 1 3* 3£, 2 s +1/2 , 3g 

17 (H u iT 2 L^)(H d )(iT 2 H u )(H u ir 2 L) 1 1* l£, 2* 1/2 , 2^ 1/2 

18 (H u ir 2 TL^(rH d )(iT 2 H u )(H u iT 2 L) 1 3* 3 L , 2 R 1/2 , 2£ 1/2 , 1* 1 L 

19 (HjT 2 JS)(H d )(iT 2 TH u )(H u iT 2 rL) 1 1* 1 L , 2 R 1/2 , 2 L 1/2 , 3* 3 L 

20 (H u ir 2 r a L-)(r a H d )(ir 2 r b H u )(H u ir 2 r b L) 1 3 R , 3fr, 2* /2 , 2 L _ 1/2 , 

21 (L^r 2 r a L)(H u ir 2 r a )(r b H d )(H u ir 2 T b H u ) 1 3^, 2; i/2 , 3^ 

22 (L c iT 2 T a L)(H d iT 2 r a )(r b H u )(H u ir 2 T b H u ) 1 3 S +1 , 2^ /2 , 3^ 

23 jL c ir 2 rL)(H u ir 2 r)(H u )(H d iT 2 H u ) 1 3 S +1 , 2 s +1/2 , lg 

24 (mr 2 T a L)(H u i7^T a )(T b H u )(H d \r 2 T b H u ) 1 3 S +1 , 2 s +1/2 , lg 

25 (H d ir 2 H u )(L-ir 2 )(rL)(H u iT 2 rH u ) 1 lg, 2^ 1/2 , 2* /2 , 3* +1 

26 (H d ir 2 r a H^(L-ir 2 r a )(r b L)(H u ir 2 r b H u ) 1 3g, 2j 1/2 , 2j 1/2 , 3 S +1 

27 {H u ir 2 L-){\r 2 H d ){rL){H u ir 2 fH u ) 1 1* 1 L , 2j 1/2 , 2j 1/2 , 3 S +1 

28 (H u ir 2 r a L^)(ir 2 r a H d )(r b L)(H u ir 2 T b H u ) 1 3« 3£, 2j 1/2 , 2£ 1/2 , 3 S +1 

29 (H u iT 2 L_ c )(L)(iT 2 rH d )(H u ir 2 rH u ) 1 1* l£, 2; i/2 , 3 S +1 

30 (ff M irVL c ) (r°L) (irV (H u ir 2 r b H u ) 1 3^, 3^, 2^_ 1/2 , 3^ 

31 {TAT 2 T a H d ){iT 2 T a H u ){r b L){H u ir 2 r b H u ) 1 3^, 3*, 2j 1/2 , 2j 1/2 , 3 S +1 

32 {L^ a H d ){r a L){ir 2 T b H u ){H u ir 2 T b H u ) 1 3^, 3*, 2; 3/2 , 3^ 

33 _{L^T 2 fH d ){ir 2 fH u ){H u ){H u iT 2 L) 1 3^, 3*, 2* 1/2 , 2j 1/2 , 1 L , 1* 

34 (L^T 2 T a H d )(ir 2 r a H u )(r b H u )(H u iT 2 r b L) 1 3^, 3«, 2*, 2* 3 L , 3* 



Table 6: Decompositions of the d = 7 operator LLH u H u H d H u at tree level. 
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and left-handed fermions, respectively, where the right-handed ones can also be represented 
by left-handed Weyl spinors after charge conjugation. All charged scalar fields must have 
an additional partner of opposite charge (not listed) to make a mass term possible in the 
superpotential. 



B Approximate diagonalization of neutral fermion mass matrix 



In our model the complete mass matrix including the flavor structure is given by 














v u Yn,u 


V U Y N ,12 

























V U Y N ,21 


V U Y N ,22 

























V u Yn, 31 


V u Yn,32 
















VvYn,11 




VvXn,31 








M N 













VvYn,12 


V u Yn,22 


VvYn,32 











M N p 



















M N 











-kiv d 


k 2 v u 
















M N p 








kiv d 


k 2 v u 





















k\v d 





—m^ 


V 



















hv u 


— 






\ 



(38) 



/ 

Using the fact, the the left-handed neutrinos are essentially massless compared to the heavy 
states we can exploit the usual seesaw formulas to obtain approximate formulas for the 
entries responsible for the mixing of the light states with the heavy states. The mass matrix 
of the heavy states is given by 



M, 



H 



/ 








M N 











\ 













M N p 












M N 











-kiv d 


k 2 v u 









M N p 








k\v d 


k 2 v u 












-hv d 


k\v d 





—m^ 




\ 








k 2 v u 


k 2 v u 


— 





/ 



(39) 



Neglecting the elements proportional to ki (i — 1, 2) this matrix is diagonalized by 



R 



H 



( 


1 





1 


V2 









1 







i 





1 













1 

V2 
















V 














1 



1 












\ 


























1 


1 




f 


^1 








) 


V2 


V2 



(40) 



The part of the mixing matrix connecting the heavy states with the light states is given by 



= mM^Rn 










( DiYnm 


D2Yn,12 


D' 2 Yn, ii 


D'\Yn,12 




DiY N2 i 


D2Yn,22 


D' 2 Yn,21 


D[Y N 22 




y DlYAr,31 


D2Yn,32 


D' 2 Yn,si 


D[Y N:3 2 



v / 2M Jv m e p 
V2M N ni£p 



v u v d (k' 2 pY Ny 



-k'Y, 



1 £ Af,12 



VuVdik^pYN^i-k'^N ,22) 

V2M N m i p 
VuVdik^pYN^i-k'^N^) 

\/2Mffm^p 
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(41) 



with 



v u Y N ,u v u Y Ntl2 
m = j v u Y N>21 v u Y N>22 ] (42) 
v u Y N ,3i v u Y N , 32 

D v u (M N m^ + 2kik 2 VdVu) ^ 

D v u (pM N m^ - 2k 1 k 2 VdVu) ^ 

D , v u (pM N m^ + 2fci^wj / 45 x 

1 " ' V2p 2 M^ 



, v u (M N m ( - 2k x k 2 v d Vu) 

2 y/2M 2 N m^ 



(46) 



k[ = k x -k 2 tan (3 (47) 
fc£ = fci + A; 2 tan /3 (48) 

Here we have the following correspondence to the couplings in Section 4, Eq. (33): 

c t = U' l5 , * = (U' l6 )* (49) 

which are the dominating ones for the lepton number violating processes. 
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